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Discovering dark matter at high energy colliders continues to be a compelling and well-motivated 
possibility. Weakly interacting massive particles are a particularly interesting class in which the 
dark matter particles interact with the standard model weak gauge bosons. Neutralinos are a 
prototypical example that arise in supersymmetric models. In the limit where all other superpartners 
are decoupled, it is known that for relic density motivated masses, the rates for neutralinos are too 
small to be discovered at the Large Hadron Collider (LHC), but that they may be large enough to 
observe at a 100 TeV. In this work we perform a careful study in the vector boson fusion channel for 
pure winos and pure higgsinos. We find that given a systematic uncertainty of 1% (5%), with 3000 
fb“^, the LHC is sensitive to winos of 240 GeV (125 GeV) and higgsinos of 125 GeV (55 GeV). A 
future 100 TeV collider would be sensitive to winos of 1.1 TeV (750 GeV) and higgsinos of 530 GeV 
(180 GeV) with a 1% (5%) uncertainty, also with 3000 fb~^. 


I. INTRODUCTION 

Among the possibilities of new physics at the Large 
Hadron Collider (LHC), observing dark matter is cer¬ 
tainly one of the most exciting. Despite the fact that 
dark matter makes up a sizable fraction of the energy 
budget of the universe, we remain in the dark regarding 
its identity. While there are many candidates, weakly 
interacting massive particles (WIMPs) are particularly 
tantalizing because of their weak-scale annihilation cross- 
section and the potential to see signals in collider exper¬ 
iments, direct detection, and indirect detection. Because 
direct and indirect detection are subject to large astro- 
physical uncertainties, producing dark matter at colliders 
seems to be especially crucial in discerning its properties. 

Many approaches to sweep through dark matter pa¬ 
rameter space have been proposed, including the use 
of effective operators [1-7] and simpified models [8-12]. 
Generically these approaches characterize the process 
PP XX where x is the dark matter particle and is ob¬ 
served as missing energy. An observation then requires a 
detectable particle to be radiated off of the initial state.^ 
These “mono-X” signatures include monojet [14, 15], 
mono-photon [16, 17], mono-Z [18, 19], mono-IT [20], 
mono-Higgs [21-23], mono-& [24-26], and mono-top [24- 
26]. 

In a UV-complete model of WIMP dark matter the 
simplest parameterization is to add an electroweak mul- 
tiplet to the standard model. The minimal supersymmet¬ 
ric standard model (MSSM) already provides the canon¬ 
ical examples of new fermion multiplets in the form of 
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^ One can also consider radiation off of the dark matter itself which 
can lead to qualitatively different signals [13]. 


neutralino dark matter: a singlet (bino), doublet (hig- 
gsino), and triplet (wino). In this work we study the 
higgsino and wino, however, the results are completely 
general and not intrinsically supersymmetric.^ 

In order to get the correct relic density, the mass for 
winos needs to be ~ 3 TeV and for higgsinos needs to be 
~ 1 TeV [27, 28]. While it is well known that non-thermal 
production can allow for different masses to saturate the 
relic density, the thermal value is a useful benchmark for 
which to aim. From the supersymmetric point of view, 
we also take the most conservative approach by assuming 
that the dark matter multiplet is only produced directly 



FIG. 1. Production cross sections where for 14 (100) TeV jets 
have pt > 50 GeV, are separated by Ag > 3.5 (4.0), and 
a dijet invariant mass > 100 GeV (500 GeV). The missing 
energy requirement is > 100 GeV (500 GeV). 


^ To leading order the pure bino has no electroweak interactions 
and overdoses the universe when produced thermally, so we do 
not consider it further. 
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rather than at the end of decay chains of heavier sparti- 
cles. In this scenario, it has been shown that the LHC 
will not be able to reach the thermal relic region [29] 
(see also [30]). Generally the most sensitive search is the 
monojet final state. Due to electroweak symmetry break¬ 
ing, quantum corrections split the masses of the charginos 
and neutralinos such that the decays of charginos can re¬ 
sult in the disappearing tracks signature which, for small 
enough splitting, can be much more sensitive than mono¬ 
jet [29]. It is still important, however, to understand and 
evaluate different search channels as their conclusions are 
complementary and can be affected by very different sys- 
tematics. 

In this vein, we study the reach for winos and higgsi- 
nos in the vector boson fusion (VBF) channel at the 14 
TeV LHC and at a proposed future 100 TeV pp collider.^ 
The importance of the vector boson fusion signature has 
long been recognized in the study of strong electroweak 
symmetry breaking [36, 37] and more recently, in dark 
matter searches [38-47]. In particular, electroweakinos 
can be produced via VBF, with a signature of forward 
jets and missing energy. We determine the optimized 
sensitivity of this search to wino and higgsino dark mat¬ 
ter, and also present results that can be applied to any 
model with new physics in VBF and missing energy. 

Note that similar work has been done in [30, 45]; this 
work is complementary (and agrees broadly with the con¬ 
clusions of [30]). In our work we perform matching for 
additional jets. This is especially important for the VBF 
signal because it suffers from large uncertainties from 
factorization scale choice [48]. Jet matching can stabi¬ 
lize uncertainties in the parton shower. Additionally we 
study the impact of future detector design on this search. 

The details of the paper are as follows. In Section HA 
we describe the analysis used to isolate the VBF signal. 
Section HB presents the details of the simulation im¬ 
plemented to compute the signal and backgrounds. The 
results for the pure wino and pure higgsino are shown in 
Section HI and for a generalized dark matter model in 
Section IV. Finally, Section V contains our conclusions. 

II. ANALYSIS AND SIMULATION 
A. Analysis Details 

Although there are many possible variations on the 
VBF signal, the general strategy is to tag two forward 
jets and then look for the new particles in the central 
region. The presence of these new particles may be in¬ 
dicated by missing energy, or by decay products like lep¬ 
tons, which can be helpful when considering compressed 
spectra. Again, we take the conservative point of view 


® See [30—35] for other 100 TeV studies relating to vector boson 
fusion. 


and assume that both charginos and neutralinos result 
solely in missing energy rather than considering the case 
where one is able to tag on some of the chargino decay 
products. 

The wino and higgsino cross-sections in vector boson 
fusion are plotted in Figure I. In practice, including addi¬ 
tional cuts improves the significance and will be described 
below. 

The largest backgrounds are Z{i'i')+iets, IF(£z/)-|-jets, 
tt, and QCD multijet. The W{^v)+iets background can 
contribute when the W decays leptonically and the lep¬ 
ton fails to be tagged because it is outside the detec¬ 
tor acceptance, not isolated, or too soft. Similarly the 
tt background is primarily from the fully leptonic decay 
where both leptons are missed. The multijet background 
arises primarily from the mismeasurement of a jet which 
mimics the jets and missing energy signal. Note that 
the vector boson -|- jets backgrounds can be produced at 
0{awal) and 0{a^) and both contribute comparably in 
our phase space. 

We employ the following analysis. The values for the 
cuts used in each signal are listed in Table 1. 

• There are required to be > njet jets, where jets are 
defined to have pT(j) > 45 GeV and pseudorapidity 

< \v{j)\- 

• The two hardest jets, ji and j 2 , are required to have 
PT{ji. 2 ) > Pr(jtag), and ji and j 2 are required to 
be opposite, i.e. ri{ji) • 77 (^ 2 ) < 0, to be separated 
in r] by more than Ar]{ji,j 2 ), separated in (j) by less 
than A(j){ji, j 2 ), and have an invariant mass larger 
than M{ji,j 2 )- 

• There must be missing transverse energy ia 

the event. This cut depends strongly on the dark 
matter mass and is optimized per mass point. 


Cut 

14 TeV 

Wino Higgsino 

100 TeV 

Wino Higgsino 

rijet 

2 

2 

2 

2 

\pU)\ 

5 

5 

7 

7 

PTijt.^) (GeV) 

45 

45 

75 

50 

AvUi,j2) 

3.75 

3.75 

4.25 

4.25 

A(t>{jl,j2) 

2 

2 

2 

3 

(TeV) 

2 

1 

10 

5 

(GeV) 

400 

- 700 

1100 

- 2500 

PT(iveto) (GeV) 

45 

45 

50 

50 

PT{e,p) (GeV) 

20 

20 

20 

20 

Pt{t) (GeV) 

30 

30 

40 

40 

vie) 

2.5 

2.5 

2.5 

2.5 

v(p) 

2.1 

2.1 

2.1 

2.1 

p{t) 

2.3 

2.3 

2.3 

2.3 


TABLE I. Cuts used for each signal for 3000 fb“^ of integrated 
luminosity. The cut is optimized for each mass point. 
Details are described in Section II A. 
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FIG. 2. Distributions of missing transverse energy for 14 TeV (left) and 100 TeV (right). All the cuts (averaged between the 
wino and higgsino values) of Table I are applied except for the missing transverse energy cut. The signal rates are multiplied 
by a factor of 100 and 10 at 14 and 100 TeV, respectively. 


• Central jets are vetoed, where a central jet, j, is de¬ 
fined to have PtU) > PrUveto) and be between the 
two forward jets rj^nin < vU) < ??max, where ? 7 min = 
niin( 77 (ji),? 7 (j 2 )) and ^max — max(r7(ji),?7(j2))- 

• Identified leptons are vetoed, where Table I de¬ 
scribes the minimum px and maximum |? 7 | require¬ 
ments on the leptons. Standard isolation criteria 
are applied and hadronic taus are tagged with an 
efficiency of 50%. 

Missing energy distributions are shown in Figure 2. 
The fx spectrum for wino and higgsino falls more slowly 
compared to the dominant backgrounds due to the pro¬ 
duction of heavy on-shell particles. Since the missing 
energy spectrum of neutralinos is directly tied to their 
masses, for each mass point the cut is separately op¬ 
timized. We optimized over all the cuts and find that, 
with the exception of $x^ ftie cuts are not strongly sen¬ 
sitive to the neutralino mass. 

QCD multijet is a potential background due to mis- 
measuring a jet resulting in seemingly large missing en¬ 
ergy. As these events tend to come from dijet events, 
the azimuthal cut between the tagging jets largely re¬ 
moves this configuration which is why we neglect this 
background. In the signal the jets recoil from the neu¬ 
tralinos and have no preference to be back-to-back.^ 

The significance for each mass point is computed as 

where S is the number of signal events, B is the num¬ 
ber of background events, and 7 _b,s is the systematic 


^ The azimuthal cut is not entirely uncorrelated with the other 
cuts. Decreasing j 2 ) tends to make the leading jets align 

more which preferentially selects events with higher 


uncertainty on the background and signal, respectively. 
The analysis in [49] employs a similar event selection 
and has a systematic uncertainty of 15%. If one over- 
optimistically scales this with luminosity, the systematics 
at £ = 3000 fb“^ would be « 1.5%. In this work ys is 
varied between 1 — 5% and 75 = 10% is fixed. While even 
this may be too hopeful, it provides a useful benchmark. 


B. Simulation Details 

The signal spectrum is generated using SUS- 
PECT2 [50] with tan /? = 20. Since we decouple all su¬ 
perpartners (except for the electroweakino multiplet un¬ 
der consideration), varying tan /3 has little effect on our 
study. Signal and background events are generated using 
MadGraph 5 v2.1.1 [51] and showered and hadronized 
using Pythia 6 v2.3.0 [52]. We use Delphes 3 v3.1.2 [53] 
as the detector simulation, using modified versions of the 
Snowmass cards [54-56] and Fastjet v3.0.6 [57] to cluster 
anti-ZcT jets with a radius of i? = 0.5 [58]. 

For the signal we generate both the 0{a^) and 
the 0{awa‘^) contributions. We do the same for the 
Z{vi^)+iets and lF(£jz)-|-jets backgrounds. We also simu¬ 
late the leptonic tt background. We neglect the hadronic 
tt backgrounds because the missing energy cut effectively 
removes them. Similarly we do not simulate the multi¬ 
jet background because the missing energy and Aiji cuts 
make it negligible. All events were matched up to one ad¬ 
ditional jet using the MLM showei-kx matching scheme. 
The events used were weighted between samples with dif¬ 
ferent generator-level missing energy cuts. 

We do not apply fc-factors to our signal or background 
events. The known A:-factors for 2^2 processes will 
differ from the fc-factors for the 2 —>■ 4 processes consid¬ 
ered in this study, and the signal cross section diminishes 
quickly as a function of increasing dark matter mass. 
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14 TeV 

95% 

1% 5% 

Scj- 

1% 5% 

Wino (GeV) 

240 

125 

135 

50 

Higgsino (GeV) 

125 

55 

70 

40 


TABLE II. Mass reach for winos and higgsinos at 14 TeV with 
3000 fb“^ of integrated luminosity. 


100 TeV 

95% 

1% 5% 

1% 

5a 

5% 

Wino (GeV) 

1100 

750 

600 

220 

Higgsino (GeV) 

530 

180 

180 

50 


TABLE III. Mass reach for winos and higgsinos at 100 TeV 
with 3000 fb“^ of integrated luminosity. 


Hence, the omission of /c-factors is estimated to be a mi¬ 
nor effect. 

In anticipation of future detector design we extended 
the coverage for jets up to \r]\ = 7 (from I 77 I = 5). This 
modification is useful to evaluate the impact of detector 
design on reach. It turns out, as will be discussed in 
Section III, changing the detector extent in rj does not 
affect the reach in this channel. 

In our studies we neglect the impact of pileup. There 
have been early studies indicating that the inclusion of 
pileup in simulations changes the efficacy of the cen¬ 
tral jet veto [59], however, new developments in pileup 
removal are expected to significantly mitigate these ef¬ 
fects [60-63]. 


III. NEUTRALINO RESULTS 

In this section, we discuss the results of our study. The 
mass reach for pure wino or higgsino at the 14 TeV LHC 
and a future 100 TeV collider with an integrated luminos¬ 
ity of 3000 fb“^ is presented in Figure 3. In all cases we 
calculate the significance using Eq. (1) and a fixed signal 
systematic uncertainty ( 75 ) of 10%. Since it is usually 
the case that S B, our results are not very sensitive 
to the particular choice of 75 . On the other hand, the 
background systematic uncertainty ( 73 ) is varied from 
1 — 5%, which gives rise to the bands in Figure 3. Ad¬ 
ditionally, we present the corresponding ranges in mass 
reach for exclusion (95%) and discovery (5cr) in Tables H 
and HI. 

As can be seen by comparing Figure 3 or Tables H and 
HI to previous studies of other search channels at 14 and 
100 TeV {e.g. monojet or disappearing track) [29, 30], 
the vector boson fusion channel is generally weaker in 
its ability to discover or exclude winos/higgsinos with 
masses on the order of a TeV. However, VBF will still 
serve as a useful and necessary complementary probe in 
future searches to confirm the consistency of a potential 


wino/higgsino discovery in monojet or disappearing track 
processes. 

As mentioned, for the 100 TeV studies, we extended 
the calorimeter coverage up to \ri\ =7, anticipating that 
future detectors may extend to |ry| ^ 6—7. In this sample, 
< 2% of jets have \r]\ > 5.5. This is in contrast with 
the vector boson fusion signal from strong electroweak 
symmetry breaking where a significant fraction of jets 
can have \r]\ > 5.5 [35]. 

The reason for the difference comes from the coupling 
of pure electroweakinos to gauge bosons. In strong elec¬ 
troweak VBF the dominant process at high energies is the 
scattering of longitudinally polarized gauge bosons, i.e. 
two goldstone bosons scattering to two goldstone bosons. 
Pure electroweakinos, on the other hand, do not couple 
directly to goldstone bosons^ so in VBF production of 
winos and higgsinos we can expect the interacting gauge 
bosons to be dominantly transversely polarized. 

When longitudinally polarized gauge bosons are emit¬ 
ted off of quarks they have similar energy, but lower pT 
than transversely polarized modes. This forces the for¬ 
ward jets to be pushed to higher pseudorapidity. With 
electroweakinos in the final state, longitudinal modes are 
not present in the goldstone limit which means the jets 
have relatively highp^ and lower pseudorapidity, \r]\ < 5. 
From this we conclude that extended calorimetry does 
not significantly impact this search. 

Finally, to evaluate the impact of integrated luminos¬ 
ity on this search we perform an extrapolation of our 
results to a dataset of 30 ab“^. Rather than generating 
an increased dataset for these, we simply multiply our ex¬ 
isting set accordingly. Compared to the third and fourth 
columns of Table I, the cuts for prUtag), PrUveto), and 
are increased to 100 GeV, 100 GeV, and 1100 — 3000 
GeV, respectively. For both wino and higgsino, M(ji, J 2 ) 
and A(p cuts of 10 TeV and 2 are used. The other cuts 
are unchanged compared to the 3000 fb“^ analysis. Since 
the missing energy cut is the most sensitive, we extrap¬ 
olate the distribution smoothly to higher values to avoid 
effects from finite statistics. The results are shown in 
Figure 4. 


IV. GENERAL RESULTS 

Our analysis is aimed at models where dark matter 
is minimally introduced as the electrically neutral com¬ 
ponent of an electroweak n-tuplet (n = 2 ,3 correspond¬ 
ing to higgsino and wino, respectively), which can then 
be produced in VBF processes via its interactions with 
gauge bosons. The signal of VBF and missing energy is 
therefore quite generic, and futhermore appears in many 
models of new physics. For instance, if there is a new 


® Electroweakinos can couple to goldstone bosons when they are 
mixed due to their higgsino-gaugino-higgs couplings. 
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FIG. 3. Wino/higgsino reach at 14 TeV (blue) and 100 TeV 
The bands sweep out varying background systematics from 1 



5%. 



V. CONCLUSIONS 


Producing and studying dark matter remains one of 
the main goals of the LHC and constitutes one of the 
primary targets for a future 100 TeV pp collider. Mean¬ 
while, the vector boson fusion channel is an important 
component of new physics searches, and here we analyze 
its relevance to dark matter. In studying various dark 
matter scenarios, supersymmetry provides a very useful 
set of examples. We presented a thorough study in the 
cases where the higgsino (electroweak doublet) and the 
wino (electroweak triplet) constitute the dark matter. 


FIG. 4. Wino/higgsino reach at 100 TeV with 30 ab“'^ of inte¬ 
grated luminosity. The bands sweep out varying background 
systematics from 1 — 5%. 


gauge singlet that is stable on collider time scales and 
possesses significant interactions with the Higgs, then 
VBF is an effective way to probe this type of model (see 
e.g. [34]). 

In light of the broad applicability of the VBF with 
missing energy channel, and the care that must be taken 
in generating the backgrounds, we provide exclusion and 
discovery contours for cross-sections vs. missing energy 
cuts. Figure 5 displays the cut efficiency times cross- 
section after all cuts from Table I have been applied (and 
fixed to the average values between wino and higgsino) 
as a function of the cut on missing energy. Given a simu¬ 
lated signal, these plots can be utilized in a simple fashion 
to obtain sensitivity for that model. To do so, one ap¬ 
plies Table Fs cuts to the signal and compares to e x cr 
to exclude/discover the signal at 14 or 100 TeV. 


In this work we analyzed the reach in the VBF plus 
missing energy channel. We found that the reach is 240 
GeV for winos and 125 GeV for higgsinos at 14 TeV. 
Going to 100 TeV, the respective sensitivity increases to 
1.1 TeV and 530 GeV. While VBF is not the discovery 
channel for electroweak dark matter, if hints of dark mat¬ 
ter were observed in a monojet search, the VBF channel 
would provide a crucial verification. This is analogous to 
the Higgs discovery in which all available channels need 
to be looked at to fully understand its properties. 

Since missing energy is a generic signature of mod¬ 
els of new physics with dark matter candidates, in Sec¬ 
tion IV we used the simulated backgrounds to set model- 
independent limits on cross-sections in vector boson fu¬ 
sion. 

Finally, we investigated the impact of extended 
calorimetry on the neutralino reach and found that it 
does not impact this search. Gompiling a list of require¬ 
ments for proposed detectors of a 100 TeV collider is im¬ 
portant: as evidenced in this study, the search for wino 
or higgsino dark matter can only be touched upon at the 
LHC and would benefit immensely from an increase to 
100 TeV. 
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VBF, 14TeV VBF, lOOTeV 



FIG. 5. Rates that can be excluded (95%) (solid) or discovered (5 (j) (dashed) at 14 TeV (left) and 100 TeV (right) for 
background systematics of 1 or 5 % as a function of the missing energy cut. All other cuts have been included, and in particular 
are chosen to be the midpoint between the optimized cuts for wino and higgsino. A signal systematic uncertainty of 10% and 
an integrated luminosity of 3000 fb“^ are assumed. 
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